Curvature of Levels and Charge Stiffness of One-Dimensional Spinless Fermions 
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Combining the Bethe Ansatz with a functional deviation expansion and using an asymptotic 
expansion of the Bethe Ansatz equations, we compute the curvature of levels D n at any filling 
for the one-dimensional lattice spinless fermion model. We use these results to study the finite 
temperature charge stiffness D(T). We find that the curvature of the levels obeys, in general, the 
relation D n = Do + 8D n , where Do is the zero-temperature charge stiffness and SD n arises from 
excitations. Away from half filling and for the low-energy (gapless) eigenstates, we find that the 
energy levels are, in general, flux dependent and, therefore, the system behaves as an ideal conductor, 
with D(T) finite. We show that if gapped excitations are included the low-energy excitations feel 
an effective flux which is different from what is usually expected. At half filling, we prove, in 
the strong interacting limit and to order l/V (V is the nearest-neighbor Coulomb interaction), that 
the energy levels are flux independent. This leads to a zero value for the curvature of levels D n and, 
as consequence, to D(T) = 0, proving an earlier conjecture of Zotos and Prelovsek. 
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I. INTRODUCTION 

The study of transport and response properties in ex- 
actly solvable one-dimensional (ID) correlated fermion 
systems, e. g. the spinless fermion, Hubbard, and t — J 
models, has-, long been an important field of theoreti- 
cal researchtrEj. This theoretical interest has been sup- 
ported by experimental studies that have found both 
qualitative and quantitative agreement between thermo- 
dynamic, spectral, and transport measurements in quasi- 
1D materials and the corresponding .properties calculated 
in ID correlated electron modelstililE 

Recently, several comparative numerical and analyt- 
ical studies have explored the differences in the trans- 
port properties between integrable and non-integrable 
ID modelsB~E3. Most of these studies have dealt with 
generalizations to finite-temperature, of Kohn's zero- 
temperature concepts and approached. In reference [ll], 
the concepts of ideal insulator and ideal conductor at fi- 
nite temperatures were introduced. These concepts refer 
to the temperature dependence of the real part of the 
optical conductivity a r (uj,T), which is given by 



a r (uj, T) = 2ttD(T)5(lo) + a reg (uj, T) , 



(1) 



where we have taken fi = e 2 = 1 , and e being the electron 
charge. The quantity D(T) is called the charge stiffness 
and characterizes the response of the system to a static 
electric field, within linear response theory. According 
to Kohn's zero-temperature criterion, the value of D(0) 
can be used to distinguish between an ideal insulator - 
D{0) = - and an ideal conductor - D(0) ^ 0. The 
quantity a reg {uj,T) is called the regular part of the con- 
ductivity and describes the absorption of light of finite 
frequency lu by the system. 



At finite temperatures, the classification of a given sys- 
tem as an ideal conductor or an ideal insulator involves 
the values of both D(T) and a (T) = a r {uj -> 0,T) as 
follows: (i) if D(T) > 0, the system behaves as an ideal 
conductor, (ii) if D(T) = and o~a(T) — the system be- 
haves as an ideal insulator, (iii) and, finally, if D(T) = 
and o~o(T) > the system behaves as a normal conduc- 
tor. 

Motivated by the results of Refs. we compute in 

this article the curvature of levels D n and study the effect 
of the temperature on the charge stiffness of the simplest 
ID lattice spinless fermion model, the Hamiltonian for 
which is 
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h = - J - c h + V J2^ - - §) . 



(2) 



where the spinless fermion operators c| , Cj obey the usual 
anti-commutation relations, means summation over 
nearest neighbors, hi = c\ci is the usual local number 



operator, ^ is the hopping integral (normally called t 
but here called j to stress the important exact cor- 
respondence between this model and the (anisotropic) 
Heisenberg S=l/2 chain), and V is the nearest-neighbor 
Coulomb repulsion. 

The., model (||) is solvable by the Bethe Ansatz 
(BA)cJ~E3 and has a metal-insulator transition at half fill- 
ing and zero temperature (in this model half filling means 
a fermion for each two lattice sites, i.e.,picrmionic den- 
sity 1/2). The metal-insulator transitions occurs when 
the first-neighbor Coulomb repulsion V is greater than 
J. This is in contrast to the standard Hubbard model, 
where the Mott-Hubbard transition occurs at half filliaa 
for any non-zero value of the on-site Coulomb repulsioncil 
(in the Hubbard model half filling means one electron per 
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lattice site, i.e., electron density 1). 

The optical properties of the model (§) can be deter- 
mined by studying the current operator j, which is given 
by 

J = ^2( c j c i+ 1 ~ c ]+i c i) ■ ( 3 ) 

3 

As in the case of the Hubbard chain, the commutator of j 
with the Hamiltonian (^|) is finite (i.e., non-zero) operator 
(i.e., not a c— number). This implies that the real part 
of the optical conductivity has finite-energy absorption, 
i.e., o~ re g(uj,T) =/= 0. The exact calculation of a r (tu,T) 
requires the full computation of the Kubo formula. The 
changes of the energy eigenvalues E n in response to an 
external flux <j) piercing the ring (formed by the ID chain 
with periodic boundaries) can be determined by solving 
Hamiltonian (Q) with the hopping integral modified by 
the usual Peierls phase factor: t — > te l ^l Na . 

The charge stiffness M{T) can be evaluated as a ther- 
modynamic quantifying using 

1 1 d 2 F 
D(T) = — >p n D n = =- 

where p n = Z~ 1 e~P En is the usual Boltzmann weight, 
Z is the partition function, 2D n = d 2 E n /d<f> 2 \<p = Q is the 
curvature of the level E n , j„ = —dE n /d(p\^ = Q is the mean 
value of the current operator in the eigenstate of energy 
E n , F is the free energy, and N a is the number of lat- 
tice sites. The summation in Eq. (Q) is not in general 
simple to perform, and it must be computed for a given 
fcrmionic density rif = Nf /N a , with Nf the total number 
of fermions. 

Many recent studies have provided crucial insight into 
various aspects of this general problem. First, using the 
Mazur inequality, Zotos and collaboratorscJ derived some 
analytical inequalities for D(T), valid for T — » oo. Their 
results imply ideal conductivity away from half filling 
both for the spinless fermion and Hubbard models but are 
inconclusive for the half-filled band case. Second, for the 
Hubbard modftL a recent analytical study by Fiuimoto 
and KawakamiEj, using the thermodynamic BAcj, con- 
firmed the exponentially activated nature of the conduc- 
tivity D{T) ch/(T) exp(-A M ff/T) (A MH is the Mott- 
Hubbard gapcS) at half filling, and ballistic transport of 
charge D(T) = Do — aT 2 (a is a positive constant), for 
low temperatures and away from half filling. These re- 
sults have been confirmed by Kirchnec. et al. combining 
both Monte Carlo and BA methods^. As a result of 
their study, these authors found numerical evidence that 
some of the conjectures of_Refs. p|-pT| may not be cor- 
rect. Third, a recent studyt£l of the "kicked t — V" model 
has focused on the transition from integrability to "er- 
godicity," in the thermodynamic limit. This model is a 



modification of the Hamiltonian (jSj) in which the inter- 
action term V is time dependent. Reference |2^ related 
the infinite number of conservation lawsE3 that charac- 
terize the integrable model (|^) to its nonergodicity, and 
argued that the deviation from the ergodic behavior led 
to the anomalous transport properties observed in Refs. 
^ [n]. In particular, Ref. ^ pointed out that nonergodic- 
ity implies an infinite transport coefficient and, therefore, 
ideal conductivity away from half filling. In the approach 
we develop in this paper, the infinite number of conserva- 
tion laws manifest themselves in the existence of only for- 
ward scattering among the exact many-body excitations 
of the system (called pseudoparticles below) at all energy 
scales. Fourth, in the regime — 2t < V < 2t (so that the 
interaction in Eq. (|^) is parametrized as V = J cos A) 
the behavior of D{T) has been studied by Narozhny_£t 
al. using exact diagonalization and finite-size scalingo. 
These authors found ideal conductivity behavior arising 
from a non-zero fraction of degenerate states in the ther- 
modynamic limit. 

It is worthwhile to remark that different authors as- 
sociate ideal conductivity with different physical mecha- 
nisms - conservation laws, nonergodicity, degeneracy of 
the many-body energy states, forward scattering (see be- 
low) - but at present the deep relation (if any) among all 
of them remains unclear. 

In this article, we focus on the regime (V > 2t) of 
the spinless fermion model in which the metal-insulator 
transition can occur and compute the needed mean values 
of the current operator in any eigenstate for the model 
(||), as well as the respective curvature of levels. In the 
regime of the metal-insulator transition we introduce for 
later convenience the parametrization V = J cosh A, with 
A > 0. 

To carry out part of our study, we use the formal- 
ism of the pseudoparticle operator algebra, which has re- 
cently been presented in detail for the Hubbard modeled. 
In brief, the BA solution (discussed in more detail be- 
low) can be shown to refer to an algebra of operators de- 
scribing "pseudoparticle" excitations. The energy cigen- 
states are characterized by the pseudoparticle momentum 
distributions, N cr( (q). Here c, 7 with 7 — 0,1,2,3,... 
are the quantum numbers that label the different pseu- 
doparticle branches, with N c ^(q) = for 7 > in the 
ground state. The BA gives the energy and other opera- 
tor mean values as functionals of the mtpientum distribu- 
tions N cn (q). Hence, as in references! 31 !! 32 !, we can employ 
a functional deviation expansion (FDE) which refers to 
the momentum deviations, SN cn (q) — N cr( (q) — N®(q), 
where N^(q) is the ground-state momentum distribu- 
tion. This leads to FDE expressions for the energy, cur- 
rent mean values, and other mean values, exactly as in 
a Fermi liquid. In addition to the central role played by 
the FDE expansions, there are other similarities between 
ID quantum liquids and 3D Fermi liquidso. These simi- 
larities have justifi.C|d-|the introduction of the concept of a 
"Landau liquid"l_in_ 3 J, which includes both the 3D Lan- 
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dau Fermi liquidsS and the ID Luttinger liquids^ " as 
specific cases. There is, however, an important difference 
between the Fermi liquid quasiparticles and the pseu- 
doparticles in ID integrable models: namely, the colli- 
sions between the quasiparticles are not, in general, of 
forward scattering type only, whereas in the case of the 
pseudoparticles in integrable models, there is only for- 
ward scattering. Moreover, the quasiparticles are not 
the true many-body excitations of the quantum liquid. 
Therefore, it is not to be expected that the transport 
properties behave in the same way for the two types of 
quantum liquids. 

The remainder of the paper is organized as follows. 
In Sec. [n] we treat the BA equations in the thermo- 
dynamic limit and compute the low- and high- energy 
excitation spectra, and the corresponding energy gaps. 
We show that the elementary excitations (the pseudopar- 
ticles) interact via forward scattering interaction (the 
/—function) only, and we give its explicit form. We 
show that the energy eigenvalues are fuctionals of the 
pseudoparticle occupancies and depend on the excitation 
spectra and on the /—functions. Using the BA equa- 
tions with a flux we obtain general expression for j n and 
D n , in the low-energy Hilbert subspace. We develop the 
FDE, which leads to simpler expressions for j n and D n , 
for the low energy eigenstates. In Sec. Ill we treat the 



BA equations in a finite-size system and derive the cur- 
vature of levels D n as function of V and N a , in the limits 
V = oo and V 3> t (1/V corrections), for all the eigen- 
states of the model. We use our results to discuss the 
ballistic transport of charge at finite temperature in this 
model. In Sec. IV we discuss our results and present the 
conclusions. 



II. BETHE ANSATZ RESULTS 



Our approach will be to use the general BA solution _ 
for the Hamiltonian (^) with a flux <f> through the ringS 
to compute the mean value of the current operator j n and 
the level curvature D n in a giyejieigenstate of energy E n 
using the formalism of Refs.E3cd. Our results will follow 
directly from the structure of the Hilbert spacefill of the 
model (||) . In the following, we discuss this structure in a 
finite-sized system and present a simple physical picture 
for this problem. 

It is well known that the model @ can be mapped 
onto the Heisenberg spin-one-half model by the Jordan- 
Wigner transformation. For a given canonical ensemble 
there are states with both S z = S and S z ^ S, where S 
is the total spin and S z is its projection in the z direc- 
tion. We shall present the results for those states that 
are characterized by S z = S, but similar results would 
be obtained for the states with S z ^ S (for example, the 
lowest energy state with S = 1 and S z = is degener- 
ate, in the thermodynamic limit, with the ground state 
characterized by S = and S z = 0) . 

The model (0) has particle-hole symmetry and, there- 
fore, without loss of generality we study the case where 
the fermion density n/ is less than 1/2. The gen- 
eral BA solution for the Hamiltonian (||) is obtained by 
combining the Jordan- Wigner transformation with the 
BA solution for the Heisenberg spin-one-half modelEZl. 
Our starting point is the set of BA equations with 
(Orbach para petrization) twisted boundary conditions, 
which arcOEj 




2 arctan f A, tan ) = "—^J. + ( 7 + + iZLEl 2 arctan [ A„ tan 

I 7 2 N a N a 'f-'n i N a V V 

7 =0p=l / \ 




(5) 



where the numbers A J C are called rapidities, the la- 
bel c is associated with charge, A 7 = coth[(7 + l)A/2], 
A p = coth(pA/2), 4> is the twist angle, and the symbol 
[77'p] is given by 

[77'p] = 1, for p = |7-7'|, 7 + 7' + 2 
[77'p] = 2, for p = | 7 - V| + 2, | 7 - V| + 4, . . . 
. . . , 7 + 7' 

[77'p] = 0, otherwise. (6) 

The energy E n of a given eigenstate is a function of the 
chosen set of numbers {/•? 7 }„ and can be written as 



E n = -Jsinh(A)^] ^ 



sinh[( 7 + l)A] 



, r ^ cosh[(l + 7 )A] - cos(A J c , 7 ) 



(7) 

The exact solution (j|) and (Q) of the model (||) shows 
that all its eigenstates, for a given canonical ensemble, 



are characterized by the quantum numbers J^ T , with 
7 = 0,1,2,3.... For a given canonical ensemble, the 
number of fermions Nt is related to the total number 
N Cjl of occupied I 3 C „ quantum numbers by 



= ^( 7 + l)7V, 



c,7 • 



(8) 



7=0 



The numbers I° c can be positive or negative integers or 
half-odd integers depending on whether N a + N c ^ are 
odd or even respectively (N a is taken to be even). The 
eigenstates of the model are characterized by the set of 
numbers N Ct j and by the configurations of the occupied 
quantum numbers V among their possible values. The 
iV c ,7 numbers give the numbers of c, 7 pseudoparticles 
and q = 2ttI 3 7 /N a is the momentum. For a given canon- 
ical ensemble the H _ are distributed in the interval 

L 5 7 
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-^2 i < p < rh i i -i n 

2 — J c,7 -S 2 ' ' * * * 5 Jv c,7 j 

Q 7 = N a - ^M(7, 7 ')^c, 7 , (9) 

7'=0 

and M(7, 7') = 2min{pf + 1, 7 ' + 1) for 7 ^ 7', and 
M( 7 , 7 ') = 27 + 1 for 7 = 7'. The 7 c .o numbers de- 
scribe the low energy excitations while the I Cl7 >o num- 
bers are related to theJ3A string solutions and describe 
the gapped excitationsEZl. 

A. Bands, gaps, and /-functions: cf> = results 

The numbers iV Ci7 and the momentum q introduced 
above have a simple physical interpretation. The index 
7 labels a set of energy bands which have their (pseudo-) 
Brillouin-zones and (pseudo-) Fermi points controlled by 
the relations (||) and (||). All these bands are separated 
from each other by energy gaps. These gaps, relative to 
the lower energy band c, 0, are given by 

A 7 = e c , 7 >o(0) - (7 + 1)^,0(^,0) , (10) 

where £ c , 7 (g) is the energy dispersion of the pseudoparti- 
cle and q Fcn = 2i:I^ x /N ai with F™ x being the largest 
quantum number in the compact symmetric distri- 
bution of these numbers around zero. For rif = 1/2, we 
have qFc,o = t/2 and 9fc, 7 >o = in the ground state. 



2tt 



/c, 7 ;c, 7 '(g,g') = Wc, 7 (g)* C:7;C:7 '('7,g') +v c .~ l >{q')<& cn >. c .~ l {q l \q) 

+ ^2 ^ ^c,7 " $c,-y ";e, 7 (j'gFc, T " , g) $c, 7 ";c,Y {jlFcg >',Q'), 
j=±l 7" 



(12) 



In equation ( |l2| ) $c,0;c,o(g, g') are pseudoparticle 
phase-shifts due to forward scattering pseudoparticlc- 
pseudoparticle collisionsc2l. In the Appendix, we give 
the integral equations obeyed by the phase shifts 

$c, 7 ;c, 7 '(g>g')- 

These /—function control the behavior of the correla- 
tion function exponents, as is the case of the conductiv- 
ity exponentsc3. The /—functions also control the en- 
ergy eigenvalues E n , which are written as functionals of 
SN C fi(q) and are of the form 

E„ -E Q = E {1) + E {2) + (higher order terms) , (13) 
with 

E^ =^ £ci7 (g)<W Ci7 ( (? ), (14) 

and 

E(2) =|EE ^c, 7 ;c, 7 '(g,g')^c, 7 (g)^ 7 '(g , ) I 
a 9,7 q',f' 

(15) 



The dispersion relations for the c, 7 pseudoparticles at 
half filling are 

sinhA K'sinhA f . 2 

e C:0 (g) = -J— J ^ \Jl - k 2 sm 2 (<j) 

e C)7 (0) = -JsinhA, (11) 

and their general equations are given in the Appendix. 
The parameter K' denotes the complete elliptic integral, 
which is a function of A and whose argument is kc3 (see 
the Appendix). Using a different method, the disper- 
sion e c *U[) was previously derived by Des Cloizeaux and 
GaudinEJ 

The ground state is characterized by both A r c . 7> o = 
and a compact symmetric distribution of the numbers 7 Ci o 
around zerccil. We represent this distribution by N® (q). 
According to Table || low-energy excitations can occur in 
the c, band for values of rif < 1/2. These excitations 
correspond to distributions of the numbers 7 Cj o different 
from the N® (q) distribution. We represent these distri- 
butions by iV c , (g) - Nlo(q) + SN cfi (q), where SN cfl (q) 
describes the deviations of the distribution iV Cj o(g) rela- 
tive to iV? (g). From this simple picture, we expect the 
system to behave in the same way as a free fermionic 
system for low temperature. 

There is, however, forward scattering among the pseu- 
doparticles, which is not present in a free fermionic sys- 
tem. This interaction is described by the /—function 
/c, 7;c , 7 ', which is given by 



and Eq is the ground state energy. Similar results to Eqs. 
(Ill)' 0)> an d (El) nave been derived for the Hubbard 
modeS^ and as long as S cn /N a -the number of excited 
pseudopartocles- is small the higher order terms are neg- 
ligible. 

The relation between the gapless c, excitations and 
the original fermionic problem is best understood in the 
Luttinger liquid paradignOCj. These excitations corre- 
spond to the density waves obtained by bosonization of 
Hamiltonian (^), which, in turn, are generated by the low 
energy particle- hole processes around the non-interacting 
Fermi surface. 

A different situation occurs if n/ = 1/2. For this 
canonical ensemble, Table | shows that the number of 
available I c .o numbers is equal to the number of c, 
pseudoparticles occupying those orbitals. As a conse- 
quence, no low-energy excitations can occur in the c, 
band. The lowest energy excitation requires that one 
N Ct i pseudoparticle is created, which costs a finite en- 
ergy. In this process two available states appear in the 



4 



c, band and therefore the remaining particles can now 
undergo particle-hole excitations within the c, band. 

For V ^> J the c, 1 excitation can be seen as an excita- 
tion to an "upper Hubbard band" and is related, in the 
original lattice, to the creation of a pair of two nearest- 
neighbor occupied sites. In this limit, the gap Ai has 
the asymptotic form Ai ~ V — 2J. This excitation 
process is the most important contribution for th|6-|Zero- 
temperature behavior of <j reg (u, 0) at half-fillingpa. In 
Fig. |l| we plot the c, band at half-filling and represent 
the lowest energy excitation to the c, 1 band (represented 
by a point). 

From this description, we see that there is a mapping 
between an insulator due to correlations in the fermionic 
picture and a band insulator in the pseudoparticle rep- 
resentation. This mapping holds, in this model, for all 
V > J and is not restricted to the "upper Hubbard band" 



strong-coupling argument given at 
ping holds for the Hubbard model 



-The same map- 



fa) 




(b) 
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contributions from the states with c, 7 > 1 excitations 
(these will play a major role only at half filling). 

We now make these considerations more quantitative. 
For the calculation of Eq. (|J), we need to compute the 
curvature of levels, D n . The eigenstates characterized 
by the q = 2irI c o/N a numbers and the corresponding 
eigenenergies E n are expressed as functionals of the dis- 
tributions N c fi(q) of these quantum numbers. We write 
N c . Q (q) = N° (q)+SN cfi (q), where N° (q) is a compact 
distribution of all I c q numbers symmetrically distributed 
around zero. If the deviation 5N c (q) is small, we can ex- 
pand the relevant quantities in 6N Cj o(q)- Using the BA 



equations with twisted boundary conditionsLrl, we next 



derive simple equations for j„ 

Following references 0il 
j n reads 



and D n as functions of 



the general expression for 



j n = J smli A > N cfi {q)- — - — — ^ , (16) 

{cosh A — cos[A c ,o(<7)J} 



where A c .o(g) is defined by taking the limit N a — > 00 
in Eq. (||) and A^ (g) is the first derivative of A^ with 
respect to <j> such that 



A* (9) 



lim 



dAi (2nli /N a ) 



d{<t>/N a ) 



(17) 



0=0 



We stress that Eq. ( |l6|) is the exact mean value of the 
current operator (|^) in any eigenstate of the Hamilto- 
nian (j|) of the low energy sector (only c, excitations 
included). That is, Eq. @ holds for any N c<0 (q). By 
methods eq uival ent to those used previously for the Hub- 
bard modeOcil and for small densities of excited c, 
pscudoparticles the FDE gives the following result for 
the mean value of the current operator (|16[) : 



= J2SN c , (q)j c ,o(q), 



(18) 



1.0 2.0 



FIG. 1. Plot of the filled c, band (corresponding to 
half-filling of the spinless fermions) and lowest energy exci- 
tation to the c, 1 band (creating a pseudoparticle in the c, 1 
band and leaving two pseudoholes in the c, band in accor- 
dance with Eq. (H)). 



B. D n for n; < 1/2 and T < u c ,og c ,o 



Due to the presence of the gaps A 7 the behavior of 
the curvature of levels D n at low temperatures and away 
from half filling is mainly controled by the c, excitations. 
Therefore in the calculation of D n below we ignore the 



with the spectrum j c ,o(q) given by 

(19) 



i=±i 



Qt'oil) = 1 + ^2 ./'I'- 'I- .I'll I , 
j'=±i 



v c ,o(l) = de a M[)/dq, and v cfi = v Ct0 (q Fct0 ). 

Equation (H3) has a simple physical meaning, since it 
is the sum of the velocity plus a dragging term due to the 
collisions among the pseudoparticleSpisimilar results have 
been derived for the Hubbard modelQci) . This picture is 



very similar to that occurring in a Fermi liquicO. The 
curvature of the levels E n can also be computed. For any 
low-energy eigenstate, its curvature D n is given by 
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D r , 



J 



sinh A 



E 



sin[A c , ( g )]A^(g) + cos[A e , ( g )][A* ( g )] 2 
{cosh A — cos[A C: o(g)]} 2 



sin 2 [A c , ( g )][A* ( g )] s 



{cosh A 



[A C) o(g)]} 3 



(20) 



where A^^g) is the second derivative of A^ with respect 
to 4> in the sense of Eq. ( p"7| ) . 

For small densities of excited c, pseudoparticles using 
the FDE allows us to express Eq. (^0|) as 

D n = D + SN cfi (q)D cfi (q) , 



D ° = E Vc >° [Qt'i(jQF,o) 



2 AT 

= rtofc.o) 2 , (21) 

Z7T 



where Do/N a is the zero-temperature charge stiffness 
and £ Ci o is the dressed charge. Identical relations have 
been obtained by other authors for other integrable 
modelse3. It is simple to show that Do = Jc,o(9f,o)/(27t), 
and similar, results have been derived for the Hubbard 
modefl'liilE3. The function D C fi(q) is an even function of 
q and is a combination of v C fi{q) 1 <& c ,0;c,o{<li q'), and their 
derivatives. We give its general form in the Appendix. 

Close to half filling, D C:0 {q) — > dv C fi(q)/dq is the cur- 
vature of the energy band and Dq — > 0, as computed 
by HaldanecZI. This, latter result is a consequence of the 
vanishing velocityCJ v Cj o(q) at q = 7r/2, as can be seen 
directly from Eq. jll|). This signals the zero temperature 
metal-insulator transition in this model. The spectrum of 
the curvature, D C fi(q), has a very simple form at close to 
half filling, when compared to its general expression given 
in the Appendix. It bears some resemblance to the re- 
sult for the independent fermion gas, but here dv Cj Q(q) / dq 
refers to the c, energy excitations of the many-body sys- 
tem. 



tion theory. In the limit V = oo the BA equations (|^) 
and (ffl) can be written as 



N a A> c % = 2nP cn - 



(7 + 1)0 



■ E Eb"^ 

7 '=0p=l 



E 

3'=1 



A^), ./ 1 V..- (22) 



and 



E,, 



JV c ,o 

- J E 

3=1 



cosAJ~ , 



(23) 



respectively. In this limit we see that the set of Eqs. 
(||) transform into a set of coupled linear algebraic sys- 
tems. This result shows that if V = oo only the A^'q" pa- 
rameters are of importance, in what concerns the energy 
eigenvalues, but, nevertheless, Eq. (22) is meaningfull for 
all A^^°. Considering only the c, excitations, Eq. (|2T 
can be solved exactly and reads 



2tt 



N a N a - N f 
2tt 



i c,0 



2?r 



N a N f N f (N a -N f ) 



l c,0 



(24) 



The picture that emerges from this solution is essentialy 
that of non- interacting particles, in an effective ring of 
size N a — N-fuB. For this case, the curvature of levels is 
simply given by 



III. ASYMPTOTIC BA RESULTS FOR 
ARBITRARY SYSTEM SIZES 



— r 

2N a ^ 

{<<>}« 



cosA c,0° 



(25) 



In the previous section we have considered the curva- 
ture of levels in the thermodynamic limit and we used 
the FDE formalism to compute j n and D n for the low 
energy states. We now consider the N a dependence of j n 
and D n . The calculation of j n and D n for all eigenstates 
for finite-size N a is a difficult problem in general, but in 
the limits V — oo and V 3> t, it can be performed explic- 
itly. Our starting point is the = oo solution which we 
will use to compute the asymptotic corrections of order 
1/V. 



A. V = oo solution 

The limit V — > oo is a well defined limit of the BA 
equations and therefore is suitable for classical perturba- 



and the zero temperature charge stiffness 
puted using Eqs. ( |24| ) in Eq. ( |25| ) givingc2 



be 



D(0) 



J sin(irN f /(N a - N f )) 
2N a wa.(n/(N a -N f )) 



(26) 



Therefore, as long as the system is not at half filling we 
expect ideal conductivity behavior at finite temperature, 
in agreement with the picture we developed in the previ- 
ous section for the thermodynamic limit. 

The set of equations (^2|) can be solved for arbitrary 
occupancy of the numbers N cn and this solution can be 
used as a starting point to an asymptotic calculation of 
the energy levels E n and its corresponding curvature of 
levels E n as function of the system size N a . The curva- 
ture of levels D n depends in a crucial way on the flux 



G 



dependence of the parameters A J r . We now determine 
this dependence by solving Eq. (|22|). 

From Eq. (||) the highest possible occupied band is 
that characterized by 7/ = Nf — 1. Let us consider an 
arbitrary occupancy of the numbers iV C:7 compatible with 
Eq. (||) . This configuration is always of the form 

N C)0 = ^0; N cA = vi , . . . , ; N c ,m = v m (27) 

and 

N cA+m = 0; . . . ; N cnf = , (28) 

with the possibility that m = 7/. In order to solve our 
algebraic coupled linear systems we first need to obtain 
A;J'°°. Incidentally we only need to compute the sum 



( 29 ) 

which, as can be seen from Eq. (p2[), is present in all 
the equations that define A^?°, with 7 < rh. Using the 
equation that defines A^'°°, we obtain 

(N a -2N f + 2(fh+ l)z/ ro ) J2 K% = 4m 
m — 1 v i 

+ Vih {l + m)0 - */ m 2 E 2( 7 + 1)A*~ . (30) 
7=0 3=1 

Substituting the above sum in Eq. ( ^2| ) for a given 7 < m 
we obtain 



(N a - 2(7 + 1 W)A# = 2< 7 + + E EI77'P] £ - 



2(7+1) 



JV„ - 2iV/ + 2(m + 1) 



Y) 27r E4- - v - E E) 2 (7+ i)Aj; 



7'=0 J = l 



(7 < to; m > 0) . 



(31) 



It is simple to prove that if only the c, excitations have 
non-zero occupancy then Eq. ( |3l| ) gives Eq. (|2^). More- 
over we see that the flux term, due to the presence of 
c, 7 > excitations is modified as follows 



4> 



e// = (Na - 2N f )cp 

N a - 2N f + 2(m + 1) 



(32) 



The above equation shows, in the regime V 3> t, that 
if c, 7 > excitations are allowed, the flux felt by the 
several c, 7 excitations is <l> e ^ , with the exception of the 
c, fh excitation. Even more important is the fact that 
the effective flux is zero at half filli ng. T his prop- 

erty turns out to be important in Section III B where we 
consider the asymptotic solution 7>t 



B. V 2> t solution: J 2 /V corrections 



Although the set of equations (|22|) has a solution for 

V = 00, the energy of the eigenstates with non-zero oc- 
cupancy of the c, 7 > bands is infinite. In the regime 

V 3> t the energy of the eigenstates with finite occupancy 
of the c, 7 > bands is given by 



-jE cqsA - 



c.O 



J cosh A y 7 AT, 



C,7 



(33) 



with A^ given by Eq. 



(31). Therefore, up the order 
where J 2 /V corrections are neg lcctcdR we have two dif- 
ferent situations: (a) away from half filling, the c, holes 
that appear in the c, band due to the creation of c, 7 > 



excitations fell the effective flux $ e /^ but they still be- 
have as independent particles and therefore the system 
display ideal (infinite) conductivity at finite temperature 
(this analysis agrees wit h th e thermodynamic limit solu- 
tion consider in Section II B); (b) at half filling & e ff = 
and because, up to this order, A J C is independent of A J C TO 
all the derivatives of Eq. ( |33"|) in order to the flux (f> are 
identically zero and the system behaves as ideal insulator 
(or at least does not display infinite conductivity see Eq. 

©)■ 

We now compute the 1/V corrections to Eqs. (|5|) 
and (Q). To keep the calculations as simple as possi- 
ble we adopt a low temperature scheme by considering 
a single tyoe of gapped c, 7 excitations, namely the c, 1 
excitationcj, which correspond to the gap A x = V — 2 J 
(negleting 1/V corrections). We show below that this ap- 
proach leads again to an ideal insulator behavior at half 
filling. We stress that our approach is equivalent to that 
of Fujimoto and KawakamiEd, for the Hubbard model, 
using the thermodynamic BA. 

Our starting point is the result for A^ and A^ 1 for 

V = 00 which we have denoted by A^'^. In the limit 

V >( we look for a solution to the parameters A J C 1 

of the form A^ 01 = A*^ + A j c ')/J (with V = J cosh A), 

where A^'g x is the 1/V correction to A^ 1 . By straight- 
forward expansion of Eq. (]?]) we obtain 



E n = - J e 



cos 



J cos 2 A{ 



cfl 



V 



•V.n ^ sm -V.U 
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_ j2 ^^K^ + j2(v _ v - 1)Nc ^ 

3=1 

The above equation tell us that we need to obtain A^^ V 
to compute E n consistently up to order 1/V. From Eq. 
(|J) wc obtain 

C,0 = ~y Sln A c,0 " ^7 2^ Sm ( A Cl - A c,0 ) • 

a 3=1 

(35) 

From the study of the A ^y solution of the previous Sec- 
tion it follows that Eqs. (|34| ) and (^5| ) only depend on the 
flux </> through the parameter A^'^°, which in the general 
case corresponds to A J C '°°. 

We now prove that the terms Y^ 1 sm ( A c'i° — A c o°) anc ^ 
J 1 C0S ( A c,i°) are indeed flux independent at half filling, 
although is not. This would mean that any ther- 

modynamic BApGalculation of the type developed for the 
Hubbard modeE3 would give a zero finite-temperature 
charge stifness. We start by obtaining A^'f 3 at half fill- 
ing. For the case we are considering this is simple to 
obtain (we give the general solution for A£°° in the Ap- 
pendix |B|) and reads 

3=1 3=1 

(36) 

Using the above result together with Eq. (^|) , it is simple 

to see that both J2? MKf " A c,o°) and £f cos(A J c ;f) 
are zero at half filling. 

Hence we see that the curvature of levels D n is identi- 
cally zero and therefore the charge stiffness D{T) is also 
zero. We remark here that our result does not dependent 
on fact that only c, and c, 1 excitations have been con- 
sidered. The main conclusion is that for arbitrary values 
of the numbers N c>ry , consistent with the Eq. (||), the 
highest excitation band c, m is always completely filled 
(see Table | ) and, in addition, is the only band that is 
flux dependent (up to order 1/V). As result we always 
have D n = and a zero value for D{T) at half filling. 



in the thermodynamic limited. Furthermore, our results 
(34^how that the Mazur inequality —D(T) > 0- derived in 
Ref. holds as an equality. 

Away from half filling, D(T) is mainly controlled by 
Dq for low T (that is, T <C v c ,o<Z-Fc,o)- Since the states 
with finite c, 7 > pseudoparticle occupancy occupied 
have energy gaps relative to the ground state, only the 
c, pseudoparticles contribute to the charge transport 
at finite but small temperatures. This implies that for 
low T the model can be thought as a gas of c, pseu- 
doparticles, and therefore D(T) is finite and the system 
is an ideal conductor in agreement with Ref.Ej. Of course, 
the curvature of levels D n is not as simple as in a true 
fcrmionic gas, due to the forward scattering among the 
pseudoparticles (these scattering processes arp-controlled 
by the phase shifts as in the Hubbard modelcJ). Never- 
theless, the forward scattering among the pseudoparticles 
will only renormalize the b coefficient in D{T) = Dq — VT 2 
away from its non-interacting value. This picture is sup- 
ported by the analysis of the BA equations (^) in the 
extreme correlated case of V — 00, which, for large N a , 
can be consider as free fermions in an effective lattice of 
size L e ff. = N a -Nf. 

Our results also provide insight into spin transport in 
the spin-1/2 Heisenberg chain, and similar conclusions to 
those above hold for the temperature dependence of the 
spin stiffness. The regime V > J corresponds to the Ising 
symmetry of the spin-1/2 anisotropic Heisenberg chain. 
In this sector, the exchange coupling in the z direction is 
larger than the coupling in the xy plane. The half-filled 
band (non-half filled band) case in the spinless fermions 
model corresponds to the zero magnetic field (finite mag- 
netization) case in the Heisenberg chain. 

As a final comment, we believe that the anomalous 
transport properties exhibited by some of the models 
solvable by the BA arise from the existence of only for- 
ward scattering among the pseudoparticles at all energy 
scales. As discussed in Sec. |, this is a consequence of 
the integrability in these models and is not expected to 
occur in non-integrable many-body systems. 
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APPENDIX A: GENERAL EQUATIONS FOR 
THE PHASE SHIFTS, ENERGY BANDS, AND 

Do,o{l) 

In this appendix we give the general equations for the 
phase shifts 4> Ci7;Ci7 /(g, q') and for the bauds e an {q). Fol- 
lowing standard methods (see, e.g., Refo) we obtain 



4>< 



- c,7;c,7 



p=i k 



arctan A„ tan 



A 



c,7 A-c,7 



^p $ c,0;c, 7 '(z, A c ,y) 



1 + - (Al - 1) cos(a; - A Cj7 ) 



(Al) 



and 



J 



dx 



sinh[(7 + 1)A] sinh A 
~ cosh[( 7 + 1)A] - cos[AO j7 (g)] 

sinh 2 \sm(x)<b C}0 ; Cn [x, A° 7 (g)] 
{cosh[(7 + 1)A] — cos(x)} 2 



(A2) 



where <J> c , 7;c , 7 '(<j, q') = $ c , 7 ;c,7'[A° 7 (g), A° Y (q')], x = 
A® (q FCt0 ) and A° (g) is A Ci0 (g) in the ground state 
and it can be written in terms of the phase shifts 
$c,7;c,7'(9,?') as 



2 arctan A~ tan 



KM 



2A 



X ° ^ C ,0;e,7'[».A° 7 ,( g )] 

-»o l + ^-(^-l)cosa: 



(A3) 



These three equations can be solved in closed form. 
In general, only a numerical solution is possible. For 
rif = 1/2 one has qpcfi = 7r /2 and A® Q (qF C ,o) — and 
this set of integral equations can be solved by Fourier 
series. The results for the phase shifts and for the bands 
are 



$CO;c,y[<o(?).A°y (?')] = 

g £ IOtV] sm[n(A o (g) _ A o ( / )}] 



-,—pn\ 



-2nX 



(A4) 



n— 1 p— 1 

and 



£ c,o(<?) 

£ C:7 (g) 



- J sinh A — J y —. — -r — 

/ — 1 cosh(nA) 

n— 1 v ' 



-J sinh A . 



(A5) 



respectively. The relation of these expressions with 
Eq.(|i"l|) axe obtained by using the delta and arnplitude 
functionsE 2 ! The parameter k is a function of AE3 k = 
4 [£„=i exp (-A(n - 0.5) 2 )] 2 / [l + 2 £ n=1 exp (-An 2 )] 1 
Both the phase shifts and the velocities enter in the gen- 
eral expression for D C fi(q) (see Eq. (f2l))). The general 
form of D C: o(q) reads 



DM = j- [vMiQtoM) 2 } +v c ,o(q)W^'°(q) 

{ 2v cA < lFcfi)^t,0-,cflU'lFc,0,q)Qt:o(j ( lFc,o) 

J=±l 

d_ 
dq~' 



q—3QFc,o 



+ J ^[v c ,o(q')^ c ,o-,cAq,q')[Q^(q')} 2 
+ v cfi (q Fcfi )$ c . ; Cfl (jq Fc , ,q)W^°(jq Fc fl) \ . (A6) 

The functions W*p°(q) and $£ 7;Ci7 ,(g, £?') are given by 



j±l 



(A7) 



9'=i«Fc,0 



and 



$ c,0; C ,o(<?:9 ) = <9c,oW ) 

X! jQtflti<lFc,o)®c,0;c,oti<lFc,0,q') 
J=±l 

d$c : 0;c,o(<?, ?') 



(A8) 



9 = J9Fc,0 



respectively. 



APPENDIX B: GENERAL SOLUTION OF A J C 
FOR ALL FILLINGS 



The equation that defines AJ°° is 
N a Al% = 2TrIl 1 + (rh+l)cj ) 



N 



7 '=0p=l j' = l 



(Bl) 



and can be cast in the form 



A 1,00 , A 2,00 



a A i,oo 



A; 



6 6 
where a, 6, and -B(0) are given by 

a = N a - 2Nf + iy rn + 2fh+l, 



b = 2m+l. 



(B2) 



(B3) 



(B4) 



and 







m— 1 ^ T 



B{4>) = (m + 1# - £ £ 2 ( 7 + X ) A Cn 



(B5) 



7=0 j'=l 



This linear system of order Vm can be solve, and the so- 
lution A;J'°° is given by 

2?r "** / 

K%W/b(a/b + Vffl -2)-v m + l) = --Y J Il,m + 
?!^(a/6 + ^-l)---^-(a/6-l). (B6) 



At half filling we have a — Vfa + 2m + 1 and A^'°° is equal 
to 



a j,oo 27r/ c,m 2TO + 1 7T ^ •/ , 5(0) 



m + 1 (^ 7l ) 2 f-< c ' m 2(m + 1)1/^ ' 

j'=l 

(B7) 



It is easy to see that the above equation gives Eq. ( pq ) 
for m = 1 
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<2f c ,0 


<o 


N c , 


dFc,7>0 


N c,~f>0 N c ,j>0 


G.S. N a - Nf 
Ex.S. 7 N a - Nf + ( 7 - 1)2V C , 7 


N a - 2N f 
N a - 2N f + 2~fN cn 


N f 

N f - ( 7 + 1)2V C , 7 


N a - 22V) 
N a - 22V/ + 2V C , 7 


N a - 2N f 
N a - 2N f 2V C , 7 



TABLE I. Available and occupied quantum numbers I c , 7 for the spinless fermion model. In first column G.S. and Ex.S. 7 
stand for the ground state and for the excited states with 2V Ci7 >o pseudoparticles in a single 7 band, respectively. In the first 
row, d c>7 gives the available J Cj7 quantum numbers, 2V^ 7 gives how many J Cj7 are unoccupied, and 2V Cj7 is the number of c, 7 
pseudoparticles. The numbers N a and Nf are the number of lattice sites and the number of fermions, respectively. 
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